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Abstract 

The spinless Salpeter equation presents a rather particular differential opera- 
tor. In this paper we rewrite this equation into integral and integro-differential 
equations. This kind of equations are well known and can be more easily han- 
dled. We also present some analytical results concerning the spinless Salpeter 
equation and the action of the square-root operator. 
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I. INTRODUCTION 



The Schrodinger equation is a very well defined partial derivative equation since its 
differential operator is a Laplacian. Moreover it reduces to a simple differential equation 
when the interaction is central. This equation has been intensively studied and it is well 
understood. A simple relativistic version of the Schrodinger equation, the Spinless Salpeter 
Equation (SSE), presents more difficulties. This equation is, 



P 2 + ml + Jp 2 + m\ ) *(f ) = (e - V(rj) *(f ), (f ) 



where mi, m2 are the masses of the particles, p is their relative momentum, V(f) is the 
potential interaction and E the eigenenergy of the stationary state ^(r ) (% = c = 1). p and 
fare conjugate variables. Actually, this last equation is not so well defined since the kinetic 
energy operator is a nonlocal one. Its action on a function /(f) is known only if /(f) is an 
eigenfunction of the operator p 2 . In this case we obtain, 



p 2 + m 2 f(r) = ya + rn? /(f), (2) 

where a is the corresponding eigenvalue of p 2 . Consequently this operator is difficult to 
handle. However the SSE is not a marginal equation. The correct description of the bound 
states of two particles is achieved with the Bethe-Salpeter equation This last reduces to 
the SSE 0, p. 297] when: 

• The elimination of any dependences on timelike variables is performed (which leads to 
the Salpeter equation ||), and 

• Any references to the spin degrees of freedom of particles are neglected as well as 
negative energy solutions. 

Despite the presence of a so particular operator, the SSE is often used (see for instance 
Refs. pHTOl) since its numerical resolution is very easy for bound states (see for example 
Refs. [p|pTHT3[| ) . However, it would be interesting to reformulate this equation in a way to 
better understand it or even to better use it. The main problem is the nonlocality hidden 
inside the kinetic energy operator. The idea is thus to extract this nonlocality to put it into 
evidence. That is why in this work we rewrite the SSE as integral or integro-differential 
equations where the nonlocality is then explicit. Moreover these kinds of equations are well 
known and defined. Another integral equation, different from those presented here, was 



found in Ref . |TT| . This formalism allows to define the action of the square-root operator on 
wave functions and calculate the resulting functions. 

This paper is organized as follows. In Sec. [H|, we present the integral and integro- 
differential equations and we calculate their kernels. In Sec. |TT|, we find some analytical 
results concerning the action of the square-root operator and solutions of the SSE. At last, 



we present, in Sec. [TV], a summary of this work. 
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II. INTEGRAL EQUATION FORMULATION 



It is possible to rewrite the SSE into different forms. The first one that we present is an 
integral equation. It is obtained using the Green functions of the kinetic energy operator. 
The corresponding kernels are different according to the particles are identical or not. The 
second one is an integro-differential equation. It is obtained taking the square of the square- 
root operator and is only valid for two identical particles. For these formulations of the SSE 



only local interactions can be used contrary to the formulation proposed in Ref. Jl 1 



A. First form 

In the introduction we gave the operator expression of the SSE. We are going now to 
formulate the integral expression of this equation. To obtain it, we use the Green functions 
of the kinetic energy operator. We first present the equal masses case. Let us consider the 
following function, 

1 r exp ( -ip.A ) 

with A = f — f and A = |A|. Then, G(A) is the Green function of the equal masses 
square-root operator since, 

2y/p 2 + m 2 G(A) = J exp (-ip.A) dp = 5 3 {f- f). (4) 

The solution of the SSE can be written as, 

9(f) = 9 (f) + J G(A) (E - V(f')) 9(f) df, (5) 
with \[p 2 + m 2 9o(f) = 0. The solutions of this last equation are 9o(f) = and 9o(f) = 



ii(mr)Yi m (r), where i\(x) = J^I l+ i(x), I v (x) is a modified Bessel function |L4|, p. 952], 
Yi m {f) is a spherical harmonic, r = \ f\ and r = f/r. In the appendix |B| we show that 9o(f) = 
ki(mr)Yi m (f ) is not a solution, the function ki(x) being defined as equal to K l+ i (x), 
with K v (x) a modified Bessel function |T^, p. 952]. Thus only the vanishing solution 
is relevant for physical problems. One sees immediately that the utilization of nonlocal 
potentials is impossible if we want to obtain a 1-dimensional integral equation. Calculating 
the 3-dimensional integral (^) we obtain 

G ^ = £k Ki{mA) - (6) 

If we consider only the case of central potentials, the wave function can be written in the 
form 9(f) = Ri(r) Yi m (f) where Ri(r) is the radial part of the wave function. If, moreover, 
we place the z axis along f, the angular dependence of A is just given by the angle between 
the vector f and the z axis. Then knowing that [15, p. 158] 
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(7) 



Eq. (H) is written as 



m 

471-2 
2?r 



y r' 2 (E - V (r')) flj (/) dr' y sin 



A 



o 



(8) 



We can now perform the integration over the angular variables. We have the well known 
relations 



2tt 



Yi m (r) df' = 2nYi m (r) 5 m0 , 



(9) 



and [15], p. 133], 



Y lm (9,tp) = exp(imip) 



21 + 1 {l-m)\ 
\\ An (1 + m)\ 



Pr (cos 6), 



(10) 



where the functions -P; m (x) are the Legendre functions Jl4, p. 998]. Finally, using Eqs. (|8|-|T0D 
and the regularized function u/(r) = rRi(r), we can write the spinless Salpeter equation on 
the form of the following integral equation: 



uAr) = — / Qi(mr, mr') (E — V(r')] «/(r') dr', 
2tt Jo v 



with 



/*7T 

Qi(mr,mr') — mrr' / sin# 

Jo 



, K i(mA) 
A 



Pi (cos 9') d9', 



(11) 



(12) 



where the functions Pi(x) are the Legendre polynomials [ III p. 1025]. The func- 
tions Gi(mr,mr') are analytic for each value of I. Performing the transformation y = 
\/r 2 + r' 2 — 2 rr' cos W, we obtain 



rr+r' / 

Qi(rnr,mr') —m K\(my)Pi 

J\r— r'\ V 



^2 j 2 

2rr' 



dy. 



(13) 



For each value of I, we can find the primitive and thus calculate the kernel. Indeed, we have 



the following relations |16|, p. 87]: 

x v+1 K v (x) dx = — x v+1 K u+ i(x), 
J x~ u+1 K v (x) dx = - x~ u+1 K u ^(x). 

With these ones it is easy to show that, 



(14a) 
(14b) 
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J x v K^ix) dx = — x v K^i(x) — (u + /i — 1) x v 1 K^(x) 

+ {v + n- 1) {v-\l- 1) J ' x u - 2 K^x)dx. (15) 
This last relation is useful when one calculates the kernel for I > 2. The modified Bessel 



functions recursion relation is [[14], p. 970] 



K, 



v+l 



X 



KJx) + K v -i(x) 



(16) 



With Eq. (p!3|) we see that the polynomial under the integral is composed only of even power 
of the integration variable. Then with relation (JT^), we can calculate analytically the kernel 
for each value of /. Another way (more simple) to calculate the kernel is to use the relation 
derived in Ref. ITTI 



Qi{x,x') 



( 1 d \ 1 

2 l z l+1 -— -[(y-z) L zK l (Vy^)-{y + z) L *K l (Vy+^) 



with y — x +x , z — 2xx'. The kernel have as expression for I = and 1 = 1: 

Go(x, x') = K (\x — x'\) - K (x + x'), 



(17) 



(18) 



Gi(x, x') = Kq{\x — x'\) + Kq(x + X 

1 



+ 



XX 1 



x + x') K\{x + x') — \x — x'\ Ki(\x — x'\) 



(19) 



The expressions (p!8|) and (|19D show that kernels decrease exponentially when arguments 
are important enough ((x + x') and \x — x'\ ^> 1) since K v (x) ~ J2/{tix) exp(— x) when 



|x| ^> 1 0, p. 961]. This behavior indicates that when the masses of particles are large, 
the nonlocality becomes negligible, the SSE becomes almost local and its spectrum coincide 
with the Schrodinger spectrum. However the explicit limit cannot be done in this formalism. 
One can remark that these kernels present a logarithmic singularity for r = r' . 

We treat now the unequal masses case. After the previous calculations, it is obvious that 
the Green function for the unequal masses kinetic energy operator is given by, 



(5(A) 



(27T) 



cxp 



p z + mf + Jp 1 + m 



(20) 



This equation becomes, if mi ^ m 2 
1 1 

GA 



p z + mf — yp z + m 



,2 > , < 2 — ' "~ ••<- \ exp ( — /'//. Aj dp. 



(21) 



(27r) 3 m\ — m\ 

The solution of the SSE can be written as a similar form than Eq. (151) and the function ^§{r 
is here always null. One can find it easily by acting the operator 



ip A + mf 



V 



ml 



on the equation Jp 2 + mf + Jp 2 + m\ I \I / o(^*) = 0. Extracting the operator (—A + m ) 
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from Eq. (pT|), performing the integration over the momenta and using the Green's theorem 
(see Ref. one can show that 



mi 



2ir 2 (ml — ml) 



J ATi(miA) + ^ ^ _ y ^ djfl + ^ ^ ^ 



(22) 



Again, placing the z axis along f and repeating the previous calculation we obtain 
1 



ui(r) 



tt (ml — m 2 .) Jo 



Qi(m\r, m\r ,s ) \E — 7(f)) Ui(r') dr' + (mi — > m 2 ), (23) 



with 



Qi(mr,mr') = Qi(mr,mr') 



d 2 1(1 + 1) 2 



dr' 2 r 
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(24) 



B. Second form 

This second rewriting of the SSE leads to an integro-differential equation. By acting the 
kinetic energy operator on the left of the equal masses SSE we obtain, 

4 (p 2 + m 2 ) *(f ) = E(E- 7(f)) *(f ) - 2 2 + m 2 7(f) *(f ). (25) 

Thus, we must calculate the action of the square-root operator on the product of functions 
7(f) \I/(f ). If we consider only central potentials this product gives a function for which the 
angular dependence is a spherical harmonic. One can show, with the formalism developed 
in Ref. |TT[| , that the radial part g(r) of the function g(r)Yi m (f), resulting from the action 
of the kinetic energy operator on the function h(r)Yi m (f), is given by 

1 f°° ~ 

g(r) = — / Gi(mr,mr') h(r') dr' , (26) 
ixr Jo 

with h(x) = xh(x). When the right-hand side integral is relevant, this rewriting gives 
meaning to the action of the kinetic energy operator on a central problem wave function. 
With this result, one obtain from Eq. (|25|) 

1 f°° ~ 

uAr) = — / Qi(mr,mr')V(r')ui(r') dr'. (27) 
2tt Jo 

This last equation is interesting because this is a kind of nonlocal Klein-Gordon Equation 
(KGE). This shows how the SSE differs from KGE when a nonvanishing interaction is 
introduced (the free SSE leads to a the free KGE with a different wave number k). Thus 
the difference is partially given by the nonlocal right-hand side of Eq. fl2"Tp. 

To show how the square-root operator can be particular, we give a simple example in 
which we show how the hidden nonlocality of the kinetic energy operator could lead to some 



dr 2 



1(1 + 1) 



m 2 + 



E[E-V(i 
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problems. Let us consider a square well with a range r c and a depth —Vq. In this case, 
according to r is less than r c or not, the right-hand side of equal masses SSE is a constant 
(E or E + Vo) times the wave function. Then one can think that in each area this equation 
leads to, 



A + Sgn(r c — r) k' 



0. 



(28) 



with 



(E + V 



E 2 



m if r < r r 



k 2 = m 2 if r > r c 

4 

The continuity conditions, for the interior and exterior solutions, at distance r = r c , fix 
the wave function and the energy. But this method to resolve the SSE does not take into 
account the nonlocality of the kinetic operator and we obtain the same kind of solution as 
the KGE. Actually, the problem is the discontinuity. The potential is not a constant for all 
values of r and the kinetic operator does not commute with it. Thus the equal masses SSE 
does not lead to Eq. (|28|). However, in Eq. fl27|) the nonlocality is contained in the potential 
part and we can now write an equation for each area (r < r c and r > r c ). We remark that 
this equation does not reduce to the free KGE when r > r c because the right-hand side does 
not vanish. Thus Eq. (|2"7f) leads to different solutions that these obtained with KGE for the 
same square well. 

We conclude this section with some remarks concerning these integral equations. The 



equation derived in Ref. |TT| is, for equal masses case 



, 2 r°° ~ 

E — V(r) uAr) = — / Qi(mr,mr') Ui(r') dr' . 
J 7T Jo 

Thus if the potential is just a constant, V(r) = V, using this equation and Eq. 
obtain, 



d 2 Kl + 1) m2 , (E - Vf 
dr 2 r 2 



4 



ui(r) = 0. 



(29) 



we 



(30) 



The nonlocality disappears and we find the corresponding Klein-Gordon equation. 

It is easy to show, with the operator expression of the SSE, that the free radial solutions 
are Ri(r) -r- ji(kr) where the function ji(kr) is the regular spherical Bessel function defined 

as 



Jiyx) 



(31) 



where the function J u (x) is the Bessel function of the first kind [14, p. 951]. The wave 
number k is given by 



k2= (i^a (mj-ml + E 2 ) 2 -ml 



(32) 
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In appendix ^ we show that, 

/ Qi (mr, mr') yr j t (yr ) dr' = yr ji (7 r ) • (33) 

Jo V7 + m 

Then one can verify that Eqs. (|TT|), fl23|) and ( p7|) are true for a vanishing potential when 7 
is replaced by /c. In appendix [B|, we show that the situation is rather different if one consider 
the irregular spherical Bessel functions. 



III. ANALYTICAL RESULTS 

In this section we present some analytical results concerning both the action of the square- 
root operator and solutions of the SSE. We first calculate the action of the kinetic energy 
operator on the functions r n exp(— mr) (n > —1 and n integer). The Fourier transform, 
f(p), of these functions are 



Hp) = (-) 



n+l 



'2 d 



n+l 



it d m n+l p 2 + m? 



(34) 



Then for each value of n we can calculate the action of the square-root operator. We have, 
for n — — 1 and n = 0, 



V 



— - exp(-mr) 2m 

A + m z = KAmr) 

7\r 



(35) 



4m 



p 2 + m 2 exp(— mr) = K (mr). 



7T 



We can get two more relation knowing the following integrals |14], p. 429] 



x 



26+1 



sin (ax 



[m 2 + x 1 



- dx = (— 



,6+1 







26+1 



2 n mi 



r(n + |) 



da 2b+l 



a n K n (ma) 



with a > 0, Rem > 0, — 1 < b < n and 



x 



26+1 



sin (ax) 



(m 2 + x 2 ) 



2\ n + 1 



dx 



\6+n 



7T 



2n\ 



0> 



dy 1 



(y b exp(-a^) 



7=m^ 



(36) 



(37) 



(38) 



with a > 0, < b < n, |arg(m 2 )| < tv. Rewriting these integrals as 3-dimensional integrals, 
extracting the square-root operator and integrating the resulting integrals one obtains the 
following formulas: 



p 2 + m 2 \ - 



gn 



dy 



- (l b exp(-r^)) 



\n+l 



2n\ 







26+1 



7cm n (2n - l)!!r dr 2b+1 



r n K n (mr) 
(39) 



and 



S 



p 2 + m z 



1 d 26+1 



(r n+1 K n+1 (mr) 



d 



_ 1 nm n+1 (2n+ 1)!! 
2n\r 



i 

- (l b exp(-r v ^) 



(40) 



These two equations are valid for < b < n. 

With these results we can find easily some analytical solutions of the SSE with a nonlocal 
interaction. Let us consider the following equation, 



2^p 2 + m 2 R {r) = E R (r) - J™ W(r, r') R (r') dr' . 



If we choose 



(41) 



W(r, r) 



a exp(— mr) 



71 



{(3 + m) Ko(mr) 



exp(—(3r') 



then using Eq. (|36D we obtain the solution 



Ro(r) -h exp(— mr), 
E n ' 



j3 + m 



(42) 



(43) 
(44) 



Thus it is easy to construct a set of solutions with the operator expression of the SSE, we 
must just consider a well chosen interaction. We can also find analytical solutions using 
Eq. ( p9|) with a nonlocal interaction: 



7T JO 

Choosing, for example, I = and 

2 

W(r,r') = Go(mr,mr') 

IT 

then we get the solution 



2 r°° ~ 

Qi(mr,mr')ui(r') dr' — Eui(r) — I W(r,r r ) ui(r') dr' 



+ 7r exp(— (3r — ar'), 



f + (m 2 - /5 2 ) 



M (r) r exp(— /?r) 
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(45) 



(46) 



(47) 
(48) 



Again we can construct with this representation of the SSE a set of analytical solutions 
corresponding to a set of well chosen interactions. 



IV. SUMMARY 

In this work our purpose was to rewrite the spinless Salpeter equation into forms more 
easy to handled, since its differential operator is so particular. To perform it we have 
extracted the hidden nonlocality of the kinetic energy operator to get it really explicit. 
Thus we have obtained in Sec. || the following equations: 
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For mi = 7T7.2, 



and, 



1 f°° / \ 

uAr) — — / QAmr^mr 1 ) (E — V (r')\ uAr') dr' ', 
2ir Jo v y 



d 2 1(1 + 1) 2 s(s-y(r) 
— m H 



dr 2 r 2 



1 r 00 ~ 

ui(r) — — / Qi(mr,mr')V(r )ui(r ) dr' 
2ir Jo 



with 



Qi(mr,mr') = Qi(mr,mr') 



d 2 



dr' 



1(1 + 1) 2 

+ - — 7T 2 - +m 2 



.12 



And for m\ ^ m 2 , 



m(r) 



7r (in 2 — m?,) Jo 



Qi(mir, m\r') [E — V(r') ) Ui(r') dr' + (mi — > 777,2). 



The kernel Qi(mr } mr') was first derived in Ref. JTT]]. It is analytic for each value of angu- 
lar momentum and is given in Sec ||. There are two integral equations and one integro- 
differential equation. These kind of equations are well known. The problem of definition, 
connected to the particular kinetic energy operator, are remove when the spinless Salpeter 
equation is rewritten into these forms or into the form (B^) (derived in Ref. [JXTJ ) . Moreover 
with this last one we can calculate the action of the square-root operator on any functions 
of the form h(r) Yi m (r) and find the radial part of the resulting functions. If g(r) is this 
radial part, we have 

1 r°° ~ 

g(r) = — / Qi(mr, mr ) h(r ) dr', 
nr Jo 

with h(x) = xh(x). When the integral in the right-hand side is relevant, this expression 
gives meaning to the action of the square-root operator on central problem wave functions. 
Indeed, this last relation allows us to explicitly calculate the resulting functions. 

In the Sec. |T1| we have found some analytical results concerning the action of the square- 
root operator on some particular functions. With these results and the integral formalism 
we have shown how to construct a set of analytical solutions of the SSE with well chosen 
interactions. 



ACKNOWLEDGMENTS 



We thank Dr C. Semay and Prof. R. Ceuleneer for their availability and useful discus- 
sions. 



10 



APPENDIX A: 



In this section we show, with the integral equation formalism, that the free radial so- 
lutions are the regular spherical Bessel functions. It is easy to prove it with the operator 
expression of the SSE because one can rewrite it as 



(A + k 2 ) *(f ) = 0, 



(Al) 



where k is given by Eq. (^). 

To prove it with the integral formulation, we use the integral expression of the kernel 
gi{mr : mr') (see Eq. (25) Ref. Q) 



m 



Qi(mr,mr') = — rr 



exp 



1 m 
u 4 



(r 2 + 



r' 2 I u 



%i ( -m 2 rr'u ) du. 



Then knowing that 

and |TJ, p. 718] 

r°° 1 / 

/ x exp(— ax 2 ) I v (f3x) JuHx) dx = — exp 
Jo 2a \ 

(with Re a > 0, Rez/ > —1), we have 



r-f\jj§i 



4a 



2a 



Qi (mr, mr') ji Oyr') dr' = 7r ji ( r yr) 

m 



3 



du 



7T 



-y/7 2 + m- 



(A2) 



(A3) 



(A4) 



(A5) 



APPENDIX B: 

We know that the relation 



V-A + m 2 G{r) = ^ ±(3 2 + m 2 G(r), (Bl) 

leads to G{r) 4- ji(flr) (plus sign) or G(r) -4 ii(fir) (minus sign). The situation is rather 
different for the irregular spherical Bessel functions n/(x) and ki{x) {ni(x) = y^iV /+ i (x), 
where N u (x) is the Bessel function of the second kind). To illustrate, we just consider the 



function k (x) = exp(—x)/x. We see that the right-hand side of the Eq. (|35| ) is different 
from zero as we could expect from Eq. (Bl|). Then we can write 



Si 



exp (-7?-) j— exp(-7r) 



/p 2 + m 2 *\ ' ' = ~ T 2 _ + f(r), (B2) 

where, 
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J 2 + m 2 + ^m 2 - 7 2 ) f(r) = 4mS(f ), (B3) 
since 

(A-7 2 ) eXp( ; 7r) = -4^(r ). (B4) 

From Eq. ( p3|) we deduce 

2 r°° p sin(pr) 

f(r) = — dp. (B5) 

vrr jo VP + m + v m — 7 

Thus we see that the singularity at the origin of the functions i%i(x) and ki(x) is very 
important and very annoying. Indeed, if f(r) was null it would be easy to solve the equal 
masses relativistic coulomb problem for 1 = 0. In this case one can show that 

p 2 + m 2 exp(— yr) = J m 2 — j 2 exp(— jr) H — , = — — — ^— (B6) 

\Jm 2 — 'j 2 r 

and, 

27 
y/m 2 — 7 



p 2 + m 2 rexp(— jr) = \jm 2 — 7 2 rexp(— 77") H — / . = exp(— jr) — 



m exp(— jr) 



i 3 

[m z — 7 2> 



(B7) 



Acting the square-root operator on the left of these relations one can verify them. The 
equation to solve is 



r 

Using Eq. (p6|), we find that the ground state is 



2 dp 2 + m 2 R^r) = (£? + -) R^r). (B8) 



i2o(r) exp ( 7T=^ r \ ( B9 ) 



with, 



^2 

And, using Eqs. (|B6|) and (]B7f) , the first excited state is given by 



V4 + « 



2m 

£ = ; (BIO) 
1 



with, 



, 16/tm \ / Km \ 
ifc r) - 1 - -t r exp - —= r Bll 



£=^W (B12) 



1+00' 



This spectrum is a rather good approximation. For example, the error is ~ 0.25%, for the 
ground state, if k = 0.456 and m = 1 GeV (see Ref. These energies were first obtained 



in Ref. [T7fl. Let us remark that the wave functions are Schrodinger-like wave functions. 
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